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Abstract. We present self-consistent solutions for the disk structure of clas- 
^ , sical Be stars. Our disk model is hydrostatically supported in the vertical di- 

rection and the radial structure is governed by viscosity (a-disks). We perform 
three-dimensional non-LTE Monte Carlo simulations to calculate simultaneously 
both the equilibrium temperature and Hydrogen level populations and to solve 
f—^ ■ self-consistently for the density structure of the disk. We discuss the general 

\^ | properties of the solution for the disk structure and test our model against ob- 

C^) . servations of S Scorpii. Our results confirm that a viscous decretion disk model 

is consistent with these observations. 

6 
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c/2 . 1. Introduction 

Observational evidence indicates that the circumstellar disks of Be stars are 
geometrically thin, optically thick disks, whose outflow velocities are small in 
comparison to their rotation speeds (see iPorter Sz Riviniusl [20031 . for a recent 
review). Dynamically this requires that the disk be pressure supported in the 
vertical direction and centrifugally supported in the radial direction. For this 
reason the disk structure is most likely to be that of a Keplerian hydrostatically 
supported disk. How material is injected at Keplerian speed into this disk is still 
unclear, but after it is injected, the gas will spread out ward (and inward if it is 
injec ted at intermediate radii), owing to viscous effects (|Lvnden-Bell Sz Pringld 
119741 V 

In a recent paper, ICarciofi h Biorkmanl (|2006al . hereafter CB) described 
a new non-local thermodynamic equilibrium (NLTE) Monte Carlo (MC) radi- 
ation transfer code that may be applied to study differentially moving three- 
dimensional circumstellar envelopes. 

As an initial application of this code, CB studied the temperature structure 
of the disks around Be stars. They found the disks to be highly non-isothermal, 
with a difference of up to 3 times between the minimum and maximum tem- 
peratures. The minimum temperatures are about 30% of the stellar effective 
temperature (T e g) and occur at the midplane typically around 3 — 5 stellar radii 
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away from the star, depending on the disk density. The maximum temperatures 
occur at the base of the star and for very dense disks can be even larger than 
T e ff due to the backwarming of the star. 

CB found that the temperature is a hybrid between that of optically thick 
disks of Young Stellar Objects and optically thin winds from Hot Stars. For the 
optically thin disk's upper layers (defined as the regions where the distance to the 
midplane is larger than one scale-height), the temperature is nearly isothermal, 
with an average temperature of about 50% of T e g. This is much cooler than the 
80% value frequently quoted in the literature (e.g.. IWaterslll986l ). 

At the midplane the temperature shows a complicated structure. It initially 
drops very quickly close to the stellar photosphere, and the temperature p rofile 
is well described by a flat blackbody reprocessing disk (|Adams et alJll987fl . 
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where T+ is the temperature of the radiation that illuminates the disk. When the 
disk becomes optically thin vertically, the temperature departs from this curve 
and rises back to the optically thin radiative equilibrium temperature, which is 
approximately constant in the winds of Hot Stars. 

The results for the disk temperature raises the question of how it affects 
the disk structure. Determining the self-consistent solution for a Keplerian disk 
is a somewhat complicated problem, because the disk temperature controls its 
geometry (via the hydrostatic equilibrium equatio n), which in turn determine s 
the heating and hence the temperature of the disk ( Kenvon Sz Hartmannlll987t ) . 



In this paper we present the first results of an updated version of the code 
described in CB. In this version of the code we self-consistently solve both the 
disk temperature and the fluid equations to obtain the disk density structure. 

In Section 2 we describe the physical properties of our disk model. In 
Section 3 we briefly describe our MC code. The properties of the solution for a 
typical Be star disk model are described in Section 4. Finally, in Section 5 we 
test our model by fitting continuum observations of 5 Scorpii. 



2. Keplerian Disks 

A pop ular model for Be star disks is the viscous decretion disk of lLee. Saio. &: Osakil 
(199!). This model is essentially the same as that employed for protostellar 
disks, the primary difference being that Be disks are outflowing, while pre- 
main-sequence disks are inflowing. In this model, it is supposed that some (as 
yet unknown) mechanism(s) injects material at the Keplerian orbital speed into 
the base of the disk. Eddy/turbulent viscosity then transports angular momen- 
tum outward from the inner boundary of the disk (note that this requires a 
continual injection of angular momentum into the base of the disk). If the ra- 
dial density gradient of the disk is steep enough, angular momentum is added to 
the individual fluid elements and they slowly move outward and form the disk. 

To critically test such decretion disk models of Be stars against observations, 
we must determine the structure of the disk, which, as stated in the introduction, 
is affected by the temperature. 
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The structure of an isotherma l viscous decreti on disk was already studied 
in the literature (see, for example, lBiorkmanlll997fl . What governs the vertical 
structure of the disk is the gas pressure, responsible for the Gaussian term in 
the formula for the disk density (written in cylindrical coordinates w, z and (j>) 

p(w,z)= e-°- 5 ^ 2 . (2) 



The disk surface density, S, is given by 

MV crit Rl /2 



S ^ = ^^lV^-lj, (3) 

This expression shows that, for an isot hermal disk, the radial stru cture is con- 
trolled by a, the viscosity parameter of IShakura &: Sunvaevl ( 1973h and Rq, the 



disk "truncation" radius. The density scale is controlled primarily by the equa- 
torial mass loss rate M, but also by the critical velocity of the star, V^tj and 
the disk temperature via the sound speed a. For large disks (i.e., large Ro), 
we have that the surface density of an isothermal disk is a power-law in radius, 
£ oc zu~ 2 . 

The disk scale-height, H, is given by 

H(w) = [a/v^w. (4) 

Since the disk is isothermal and the <j> component of the velocity is Keplerian 
(^V> = V^ritro -1 / 2 ), we find the familiar result that an isothermal disk flares as 
H oc w 1 ' 5 . Hence, it follows from Eq. |2]that for an isothermal disk the density 
falls very sharply with radius as p oc tu~ 3 ' 5 . 

To conclude this brief description of the properties of isothermal Keplerian 
viscous disks, we must describe their velocity structure. The cj) component is 
Keplerian, as stated above, but those disks must have an outflow velocity. This 
velocity is given by the mass conservation relation 

M 

(5) 



As mentioned in the introduction, CB found that the disks of classical Be 
stars are highly non-isothermal. A temperature structure affects the density 
structure of the disk via the sound speed in Eqs.|21and0J It is straightforward 
to show that the formulae for non-isothermal viscous disks are the same as 
Eqs. <PJ) to ©, with the temperature replaced by the weighted average of T in 
the vertical direction ( see ICarciofi $z Biorkman1l200 6b. for details) 

m _ S-ooT{w,z)p(w,z)dz 



3. Monte Carlo Code 



For the results shown in sections 4 and 5 we employ an updated version of the 
code of CB. The basic steps of the code are as follows: 
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1. The code starts with a constant temperature for the disk, an initial guess 
for the Hydrogen level populations and the isothermal viscous decretion 
density [Eqs. © to ©]; 

2. A full NLTE MC radiation transfer, using the Sobolev approximation for 
the line opacities, is performed. 

3. The heating, ionization and excitation rates sampled during step 2 are used 
to solve the rate equations and the energy conservation equation, to obtain 
news values for the disk temperature and Hydrogen level populations. 

4. The new disk temperature is used for calculating the new density structure. 

5. Steps 2 — 4 are repeated until all the state variables and the disk temper- 
ature have converged to a solution. 

6. If wanted, a post-processing code is run to calculate observables (SED, line 
profiles, images, etc.) using the results of step 5. 



4. Properties of the Solution 



Table 1. Parameters of the disk and the central star 



T cS (K) R*(R Q ) M(M /year) a V cii t (km/s) R (R*) Po (g cm" 3 ) 
19,000 5.6 5xl0~ 12 0.1 560 1000 3.7 x 10 11 




Figure 1. Temperature distribution. Each panel shows the temperature as 
a function of x and z for a different radial scale. 

In tabled we list the parameters of the disk model we choose as an example 
to illustrate the properties of our solution for the disk structure. The model is 
similar to the ones used by CB and consists of a B3IV star surrounded by a 
steady state disk, sustained by an equatorial mass loss rate of 5x 10~ 12 Mq /year. 
We also carried out calculations for other spectral types and mass loss rates. 
The details of the solution varies with the choice of parameters, but the general 
properties of the disk structure are similar. 
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Figure 2. Radial temperature structure. Left: Temperature along the mid- 
plane. Right: Weighted-average of T in the vertical direction [Eq. ©]. 



CB studied the temperature properties of disks with a known density func- 
tion, corresponding to the isothermal viscous disk of Eqs. © to (J5J). The tem- 
perature structure of our self-consistent Keplerian disk is shown in Figures ^ 
and|21 Interestingly, the temperature structure is nearly identical to the results 
of CB (see, for instance, their Figures 6 and 7). The disk is isothermal in the 
upper layers and also in the vertically optically thin portions of the midplane 
{w/R* > 10). The inner disk (zu/R* < 10) behaves initially as a flat reprocess- 
ing disk [eq. (Q) > thin line of Figure |2J-a] and departs form this curve when it 
starts to become optically thin vertically, which happens for w / R± ~ 3 (see CB 
for further details). 
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Figure 3. Disk density structure. Left: Opening angle. The thick line 
corresponds to the self-consistent solution for the disk opening angle [eq. Q]. 
The dashed lines correspond to opening angles for scale-heights parameterized 
as H(vj) oc vj 13 . From bottom to top, j3 — 1.0 to 1.5, in steps of 0.1. Right: 
Density profile. We compare the density calculated from our MC simulations 
(thick line) with the p oc tu~ 3 ' 5 isothermal density (thin line). 



In Figure |31 we analyze the properties of our solution for the disk density to 
understand how our self-consistent solution differs from the isothermal solution. 
In the left panel we show the disk opening angle vs. radius, defined as 
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Figure 4. Best-fitting model. Left: Our best-fitting SED (solid line) is 
shown along with the observed active-phase SED (filled circles) and the pre- 
active phase SED (stars). The other lines correspond to the scattered, emit- 
ted and attenuated stellar fluxes (dotted, dashed, and dash-dot lines, respec- 
tively). Right: The best-fitting model polarization (thick lines) and the ob- 
served polarization (thin line) are shown. The solid thick line correspond to 
the model polarization with the depolarizing effect of the secondary included, 
and the dotted thick line correspond to the uncorrected polarization. 



For reference, we also show the opening angle corresponding to disk scale-heights 
given by H(w) oc w 13 . A (3 of 1 corresponds to a non-flared disk (wedge-shaped) 
and a (5 of 1.5 to an isothermal disk. 

The opening angle for the inner disk {wjR* < 10) is very small, in agree- 
ment with the general belief that the disks of Be stars are geometrically very 
thin. However, the behavior of the opening angle is rather peculiar; we observe a 
decrease of the opening angle (de-flaring) between 1 and 2 stellar radii. This de- 
flaring is a consequence of the very sharp fall of the temperature (see Figure |2J) • 
The disk flares significantly only for m/R* > 10, where it becomes isothermal. 

The radial structure of the disk density is shown in the left panel of Fig- 
ure 01 We also show for comparison the p oc ro~ 3,5 curve corresponding to the 
isothermal solution for the density. We see that the slope of the MC density 
profile for the inner disk departs significantly from the -3.5 value; it is much 
shallower close to the star and much steeper for 4 < rn/R* < 10. This is, again, 
a result of the temperature structure in the midplane. From the observational 
point of view this is an interesting result, because the slope of the density profil e 
can, in principe, be mapped via a detailed study of the IR excess (Waters 1986). 



5. The Case of 5 Scorpii 



Table 2. Stellar Parameters for 6 Scorpii 

i?* (Rq) TeffK M* (M ) V^t (km/s) 

7 27,000 14 620 38 ±5 



To critically test our model against observations, we model continuum ob- 
servations of 5 Scorpii. This star has been a standard of spectral classification 
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Table 3. Best- fit Disk Parameters for 6 Scorpii 

M (M /year) Po (g/cm 3 ) R d (R+) i (°) 
2.5 x 1CT 9 5.2 x 10" 10 5 33 



(BO IV) since long ago, but in the summer of 2000 it started to show signs of cir- 
cumstellar activity that have lasted to present day. The object is a non-eclipsing 
binary system with a ~1.5 mag optically fainter secondary companion with an 



orbital period of 10.6 years (|Beddindll993l ) and a highly eccentric orbit (e=0.94, 
iMiroshnichenko et al.ll200ll ). We use in our modeling data corresponding to the 
SED before and after the Be ph ase. We also model t he continuum polarization 

between 3000 and 10000 A fsee ICarciofi et al.ll2006cl . for further details). 

The parameters of the primary companion were extracted from lMiroshnichenko et all 

(2001) and are summarized in Table |2J We assume that the disk around the 



primary companion is in the orbital plane of the binary system, so the disk 
inclination angle is the same as the orbital inclination angle i. 

Figure |I] shows our best-fitting SED and polarization along with the ob- 
served data. The parameters of our best-fitting model are listed in Table |3 Our 
model reproduces the flux observations quite well. A particularly significant re- 
sult is the fact that the change in slope between the visible and the NIR fluxes 
is well described by the model. Also, the model reproduces well the polarization 
in the Paschen continuum (le ft curve of Figure . 

In ICarciofi etaD (|2006c) we discuss the above results in detail. There are 



some shortcomings in our modeling that we do not discuss here, given the limited 
scope of this paper. However, the fact that our models reproduce well the 
observations for 5 Scorpii is an important evidence that this star (and probably 
the other Be stars) is surrounded by a non-isothermal viscous Keplerian disk. 
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